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Abstract. The huge worldwide investment in CMB experiments should 
make the Sunyaev-Zeldovich (SZ) effect a key probe of the cosmic web 
in the near future. For the promise to be realized, substantial develop¬ 
ment of simulation and analysis tools to relate observation to theory is 
needed. The high nonlinearity and dissipative/feedback gas physics lead 
to highly non-Gaussian patterns that are much more difficult to analyze 
than Gaussian primary anisotropies for which the procedures are rea¬ 
sonably well developed. Historical forecasts for what CMB experiments 
might see used semi-analytic tools, including large scale map construc¬ 
tions, with localized and simplified pressure structures distributed on a 
point process of (clustered) sources. Hydro studies beyond individual 
cluster/supercluster systems were inadequate, but now large-volume sim¬ 
ulations with high resolution are beginning to shift the balance. We illus¬ 
trate this by applying “Gasoline” (parallelized TreeJ-SPH) computations 
to construct SZ maps and derive statistical measures. We believe rapid 
Monte Carlo simulations using parameterized templates centered on point 
processes informed by optical and other means on the observational side, 
and by hydro simulations on the theory side, should play an important 
role in pipelines to analyze the new SZ field data. We show that localized 
sources should dominate upcoming SZ experiments, identify sources in 
the maps under filtering and noise levels expected for these experiments, 
use the RCS photometric optical survey as an example of redshift local¬ 
ization, and discuss whether cosmic web patterns such as superclusters 
can be enhanced when such extra source information is supplied. 
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1. Semi-analytic cf. Hydro Approaches to SZ Forecasts 

1.1 The Resolution of Upcoming SZ Experiments: The Compton up- 
scattering of CMB photons by hot inhomogeneous (nonlinear) gas leads to sec¬ 
ondary CMB anisotropies AT/T{q,i') = —2yc{q)'4’K{hv/ksT^) in direction q 
at frequency v, where yc = (crr/meC^) / UsksiTe — T^)dl.o.s is the Compton 
y-parameter and 'ipK depends only upon frequency.|^ It is because yc is a direct 
probe of the line integral of the electron pressure in the hot intergalactic medium 
that the heavy investment in SZ experiments is so worthwhile. 

SZ observations of individual rich clusters have been possible for a decade, 
are now routine (e.g., Carlstrom et al., 1999), and complementary to X-ray, op¬ 
tical and weak leasing observations. For resolved sources, the surface brightness 
of an SZ source is independent of its redshift. Even with the expected source 
evolution, this property should make the SZ effect a valuable probe of the clus¬ 
ter/group near/mid-field even at redshifts z ~ 1, when we expect the system 
to be in a very active merging state. Thus, the era of blank field (or ambient) 
SZ surveys is upon us, some targeting the ~ V resolution well-matched to the 
cluster/group system at z ~ 1, others with ~ 5' resolution, probing larger sky 
fractions, albeit with considerably enhanced source confusion. Specifications of 
a sample of upcoming/proposed SZ experiments are: 

Bolometer-based: ACBAR: Viper telescope, 16 element, multi-frequency, ~ 4' 
/w/im resolution, now; Bolocam/CSO: 10.4m CSO telescope, 151 pixels, 150, 220, 
270 GHz, 1', fall 2001; Bolocam on the LMT?; ACT: 3 32x32 pixel bolometer 
array, 1.7', proposed; Planck: ~ 7' at 150 GHz, ~ 5' at 220 GHz, full sky, 2007. 
HEMT-based Interferometers: OVRO mm array: 6 dish, 10.4m, 30 GHz, 
now; BIMA: 10 dish, 6.1m, 30 GHz, now; CARMA: OVRO-I-BIMA; CBI: 13 
dish, 0.9m, 30 GHz, ~ 4'; SZA: 6 dish, 3.5m, 30 GHz (-1-90 GHz), ~ 2'; AMIBA: 
19 dish, 1.2m-|-0.3m, 90 GHz, ~ 2'; AMI: 10-1-8 dish, 3.7m-|-13m, 15 GHz, ~ 2'. 
The resolution can be improved by spreading the dishes to longer baselines. 

1.2 Historical Semi-analytic SZ Forecasting: SZ estimates and limits de¬ 
rived from them were influential ever since Sunyaev and Zeldovich proposed the 
effect. First, attention was paid to baryon-dominated (BDM) models, with en¬ 
tropy injection via shocks or radiation, then to shock-heated neutrino-dominated 
(HDM) models, and then the many variants of cold dark matter (CDM) mod¬ 
els, often with strong “feedback” of energy into the pregalactic or intergalactic 
medium. For example, in the eighties one of us (Bond et al., 1980s) did HDM 
forecasts, first using the popular pancake treatment of structure formation, then 
a better cluster-based treatment using density peaks of various masses, includ¬ 
ing both Poisson and continuous clustering contributions of these “shots” (B88). 
Explosion-dominated models of structure formation in CDM, BDM models, from 
very massive objects, galaxies, superconducting strings, extreme preheating, etc. 
were also addressed in B88. All of these models became severely challenged by 


^Here, ctt is Thompson cross section, rrie is the electron mass, Ty is the CMB photon tempera¬ 
ture, and dl.o.s is the g-line-of-sight radial distance element. = 2 — (a:/2)(e“^-|-l)/(e^ — 1) 

is 1 in the Rayleigh-Jeans region, zero at 217 GHz, negative above. 
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the COBE/FIRAS data which gave a 10“^ 95% CL upper bound to the allowed 
fraction of the CMB radiation in a Compton cooling distortion. 

Early SZ maps using cluster/group-scale peaks in CDM models were also 
constructed in the eighties (B90). In the nineties, the calculations of SZ maps 
became more sophisticated with the peak-patch technique, which included cor¬ 
rect spatial clustering of halos, but “painted on” simple parameterized pressure- 
profiles within the halos (Bond et ah, 1990s, Bond & Myers 1996, Bond & 
Crittenden 2001 [BCOl]). 

1.3 Hydro Approaches: Early single-cluster SPH calculations (Bond et ah, 
1990s) became more sophisticated as the codes and computing power improved, 
as described for example in the ITP “adiabatic” cluster comparison test (Frenk 
et ah, 2000), for which SZ, X and weak leasing maps are shown in BCOl. Ideas of 
the cosmic web interconnections of clusters of peak-patches were used to create 
optimally-designed rare supercluster (TreePM-SPH) simulations (Bond et ah, 
1998, reviewed in BCOl) that included a proper tidal field acting on its 104 Mpc 
high resolution patch (comoving lattice spacing 1.0 Mpc, best z=0 resolution 
20 kpc) and cooling (but no feedback) to see which, if any, of the X, SZ or weak 
leasing probes could be sensitive to mid-field and far-held structures around 
clusters and groups (e.g., the far-held hlamentary bridges). As expected, we 
found weak leasing probes the far-held better than SZ, which does better than 
X. However, a critical issue is how to reveal such extended-source patterns, given 
the projected contributions from the clusters and environs behind and ahead of 
the supercluster targets (Sec. 3). 

In recent years, simulated SZ maps generated from hydrodynamical cosmo¬ 
logical simulations have been used to predict what the experiments should see, 
mostly with emphasis on low order statistics such as the angular power spectrum. 
Even with the great improvements in computing power and codes we have seen, 
only a relatively small number of hydro realizations per SZ map-making exercise 
are being done, and so all such work remains statistically incorrect. Nonethe¬ 
less, we believe these approximate treatments are useful steps along the path to 
that brave day when the full redshift range relevant to the projected maps from, 
say, 0 to 2, is do-able, the space tiled by contiguous simulation patches self- 
consistently constructed to have coherent long-waves joining them, all at the 
required resolution. This is routinely and rapidly done with the painted-halo 
peak-patch approach (Bond & Myers 1996). 


2. Approximate SZ Map-Making using High Resolution Hydro 

2.1 Hydro Simulations with Gasoline: We are applying the very efficient 
parallelized “Gasoline” code, developed by one of us (JW) in collaboration with 
Joachim Stadel and Tom Quinn, to a new round of huge LSS simulations target¬ 
ing the SZ effect. This treeJ-SPH code uses the smooth particle hydrodynamics 
method with a pure tree-based gravity solver and has spatially-adaptive time 
stepping. It has been parallelized using the MPI architecture and shows excel¬ 
lent scalability. The results presented here are based on the analysis of a 200 Mpc 
high resolution simulation of a “now standard” ACDM cosmology (Ha = 0.7, 
Qm = 0.3, Dbh? = 0.020, h = 0.7, rig = 1, cj 8=0.90) using a “workhorse” level of 
256^ dark matter and 256^ gas particles (lattice spacing 0.78 Mpc, best resolu- 
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tion 15 kpc). We are now analyzing a larger (400 Mpc box) ACDM simulation 
using 512^ dark matter and 512^ gas particles (now with = 0.022), per¬ 

formed on a large-memory 114 (667 MHz) processor COMPAQ SC cluster at 
McMaster, that required about 80 GB of memory and took ~ 40 days of wall 
time to run — a “state-of-the-art” simulation for SZ studies. The 256^ run re¬ 
quired 15 GB of memory and less than 6 days on 88 processors. The specific 
calculations were “adiabatic”, in the sense that only shocks could inject entropy 
into the medium. Other sources of entropy injection, from galaxy winds for ex¬ 
ample, are expected and generally may be highly inhomogeneous, a complication 
important to explore but which we ignore for this exercise. 

2.2 SZ Maps with Experimental Beams and Noise: Given only one single 
medium-sized simulation volume whose output is densely sampled in redshift, 
we create a “pseudo-realization” of the cosmic structure by stacking randomly 
translated and rotated copies of the same periodic box, an approach taken by 
other authors (e.g. da Silva et ah, 2000; Springel et ah, 2001). The SZ maps are 
generated by projecting the gas pressure in each copy of the box along the line 
of sight and later adding all the projections. This is inadequate because objects 
at various times in their history are seen again and again, but one hopes that 
it is good enough for statistical exploration if the size of the volume is not too 
small. As computing power improves, the number of boxes to be brought into 
the mix will be larger and larger, as will the size of the individual ones. 

Fig. I shows the effect of smoothing for beams spanning the range of the 
proposed SZ experiments, and also how the addition of varying levels of ho¬ 
mogeneous noise obscure the SZ source structures. The maps were made at 
Rayleigh-Jeans wavelengths, so AT/T ~ —‘^yc- Each experiment will have its 
own special Gspace filtering which truncates at low i as well as high, but for our 
purposes this simple Gaussian-filtered map will do to make our points. Note the 
high source confusion expected for experiments with beam sizes ^ 5'. 

2.3 Finding SZ Sonrces Using SZ Maps Alone: The preliminary approach 
to source identification shown in Fig. used an algorithm borrowed from optical 
astronomy developed for galaxy identification by Yee (1991), but more optimal 
object identification algorithms appropriate for separating the SZ signal from 
primary CMB signal and Galactic foregrounds as well as noise (following the “op¬ 
timal” algorithms described in BCOl) are needed to analyze scanned maps, both 
single-dish bolometer-style and interferometer-style using mosaicing/drifting. 

2.4 Where Does Most of the Contribution Come From? It is clear 
visually that most of the contribution to the noise-free SZ maps comes from 
concentrated sources. To address this quantitatively, we applied overdensity 
cuts 5g^cut to the gas particles in our SPH simulations that were allowed to con¬ 
tribute to the SZ effect, and compared the resultant maps visually and using 
the one-point distribution function of the pixel values of the temperature fluctu¬ 
ations (Fig. Q). Except for the density cut, all other aspects of the “stacking of 
boxes” map-making method were treated in the way outlined above to facilitate 
comparison. Since SZ experiments typically filter long wavelengths (some much 
more than others), subtracting off the means AT of the maps from the pixel val¬ 
ues ATp is a good way to compare what the effect of the cuts is. The histogram 
sequence with varying 5g^cut in Fig. shows that for yc above yT convergence 
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Figure 1. The left panels show a typical 2° x 2° SZ map generated 
with the pseudo-realization method of Sec. 2.2 applied to a “workhorse” 
200 Mpc Tree-|-SPH ACDM simulation, with 256^ N-body plus 256^ 
SPH particles, rather than the 400 Mpc 512^ -|- 512^ simulation. The 
maps are smoothed with fwhm filters of 1.5', near-optimal for the high 
redshift cluster/group system, and the source-confused 5' appropriate 
to ACBAR, Planck, etc. The right panels show the same 1.5' map 
after an increasing amount of instrumental noise has been applied - 
mimicking proposed AMIBA specifications and integration times from 
“deep” (300 hrs) to “shallow” (3 hrs). The 30 hr AMIBA noise level 
roughly corresponds to ~ 200 hrs for current Bolocam specifications. 
Circled objects were found in the maps using an efficient algorithm 
developed to identify optical galaxies in noisy CCD images (Yee 1991). 
If CMB primary and Galactic foreground signals are added, better de- 
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2° X 2° maps - 0.1 < redshift < 5.0 



Figure 2. Histograms of the fraction of pixels in a series of SZ maps 
corresponding to density cuts 1 + Sg^cut £ {0,20,50,100}. The sub¬ 
tracted means 2yc are {5.80,4.91,4.15,3.46} x 10“® respectively. Con¬ 
vergence beyond the means to a nearly common histogram suggests 
that a localized pressure-source model should capture much of the ob¬ 
servable range. How many of the far-held pixels (where the histograms 
deviate) should be added to adequately extend the sources depends 
upon the noise and hltering characteristics of each experiment. 


occurs for 6g^cut ^ 100. Our dg^cut SZ maps show that the pixels just below this 
convergence point, where the deviations start to appear, are contiguous to the 
high pixels, and dehne the mid-held/far-held skirt that would dehne the pressure 
prohle for the extended-source model. How far one needs to go is very dependent 
upon the noise and hlter characteristics of the various experiments that we are 
exploring. At lower levels, the pixels break out into the distributed filamentary 
gas (i.e., outside of groups.) SZ maps of thin slabs in redshift (with density 
cuts) are proof of this extended-object-contribution picture. Low mass groups 
(random superpositions) populate the interfilamentary regions in deep SZ maps, 
although they would be filtered out by experiments and not be discernible from 
the interfilamentary gas . 

2.5 Thermal SZ Power Spectra: Fig. contrasts the current CMB data on 
primary anisotropies with the SZ Ci from our ACDM SZ maps, both the small 
ones using the “stacking of boxes” and the larger scale ones using peak-patches 
(BCOl). Even though power is concentrated in sources, the Q spectrum derived 
from our SZ maps does show the angular scales we would most like to probe. In 
Fig. ^ we can see that at Planck resolution, ig ~ 1600, z < 0.5 sources are the 
most prominant contributors, but at the 1^2' scale of ground-based bolometer 
arrays and interferometers, the z > 0.5 cluster system dominates. At high i, the 
fate of gas in small clusters and groups is an ongoing uncertainty, reflected in the 
location of the Q peak and the i behaviour beyond. The overall Q magnitude 
is quite sensitive to cosmological parameters, in particular ug and its evolution. 


3. Simulating Source Localization and Removal in SZ Maps 

3.1 The Needle in a Haystack Problem: To address whether physically- 
connected supercluster concentrations can be identified in the pressure-projected 
SZ maps, we generated individual SZ maps for the largest superclusters found 
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at intermediate redshifts in the simulated box (e.g., the largest at z=0.82), and 
superimposed them on top of our pseudo-realization SZ maps. We found that 
the strong, highly concentrated, sources still stand out, but not the filamentary 
structures: in fact, many apparent filamentary sources that do stand out are 
actually due to random superposition. 

3.2 Templates from non-CMB Surveys to Localize in Redshift: As well 
as “cleaning” future SZ maps of other CMB signals and noise, we hope it will 
be possible to remove the contribution from foreground and background concen¬ 
trated SZ sources to allow access to structures in a target redshift range. This 
will not be possible without adding extra information, e.g., from deep surveys 
in the optical, X-ray and weak-lensing. The mass, redshift, size, orientation 
and other information on optically-detected clusters can be used, e.g., to place 
parameterized intra-cluster gas profiles at their locations. Although this would 
only allow a (rather dirty) form of cleaning to enhance SZ structures, in con¬ 
junction with Monte Carlo error estimations (and follow-up observations), with 
appropriate algorithm development it would be a dramatic improvement over 
what can be done with field SZ maps alone. 

We are studying the viability of this approach by coupling our simulated 
SZ maps to the cluster detection characteristics of the very successful 100 deg^ 
(22 patches X 5 deg^) optical photometric “Red-Sequence Cluster Survey” just 
completed by Yee and Gladders (2001). Proposed optical surveys such as RCS2 
or VISTA, covering 1 and 2 orders of magnitude larger areas, are planned for 
the next five years. 

When a pseudo-realization (pyramid) of the cosmic structure is generated 
by stacking copies of the simulation box, we not only create the corresponding 
simulated SZ maps, but also generate a catalogue of the clusters/groups in the 
realization by properly translating, rotating, etc. the halos identified in the box 
with a group-finder at each redshift. Fig. ^ shows how the clusters in one such 
pseudo-realization are related to the corresponding SZ map and a visualization of 
the detectability and localization of the clusters if we had a photometric survey 
with the RCS characteristics. 

3.3 SZ Maps Using Cluster Catalogues: Since the dominant contribution 
to SZ maps is due to concentrated sources, as shown above, one can also cre¬ 
ate fast Monte Carlo simulations by using catalogues of halos identified in very 
large volume (but low resolution) N-body only simulations or generated by the 
“peak-patch” technique (Bond &: Myers 1996). The intracluster gas is mod¬ 
eled by “painting” simple parameterized profiles on the halos, though with the 
improvement in progress of using (variously averaged) pressure profiles of the 
clusters identified in our high resolution boxes. These approaches offer the ad¬ 
vantage of allowing a much faster generation of cosmological realizations of the 
structures most relevant for the study of the SZ effect. Moreover, rare events in 
the cosmic structure, missing in periodic box hydro simulations, are now accu¬ 
rately included because coherent box-to-box long waves are. The drawback is, 
of course, that much hinges on profile choice. 

3.4 Summary and Future Explorations: We described here our use of high 
resolution hydro simulations to develop the “extended source” model of the 
highly non-Gaussian SZ distribution. We plan to generate more and higher res- 


Bond, Ruetalo, Wadsley Sz Gladders 


10 


102 


103 


104 


103 



0.01 - 


10 - 


T—I I M lllj-1—I M Mllj-1—I M llllj-1—I I I llll|-1—I I I llllj-TTI 

Aug. 2001 CMB data 
cf. a Boomerang best fit aka ACDM 


.[(AT/T)/10-5]2 


hydro, z = 0.1-2'(0.1-.5 .5-2) 



5' fwhm 1.5' fw\hm 
multipole L , slz slz , 

I I I I iiiir III .I I I I I Mill I . I I I I 


10 102 103 104 


103 


Figure 3. SZ power spectrum Q for one of our 4 deg^ SZ maps gener¬ 
ated from our 200 Mpc ACDM simulation. It is compared to a (scaled) 
Ci for a 25 deg^ peak-patch simulation using a slightly different ACDM 
cosmology. At 5' resolution, 2 : < 0.5 sources dominate, but at l'-2' 
resolution, the high 2 ; cluster/group system dominates. Sample vari¬ 
ance affects the hydrodynamic result at lower i, and at high i, short 
distance structure is treated differently, reflecting computational and 
feedback/cooling physical issues still unresolved in the subject. Just as 
in Fig. 1^, the range included all 2 : < 2 contributions to the TreeJ-SPH 
computation, but the adiabatic calculation ensures that big galaxies 
(at the resolution edge) and very small groups contribute, especially at 
higher i, probably more than they should. The peak-patch groups had 
a mass cut imposed, a sharp-threshold model of group winds motivated 
by this uncertain physics: hence the smaller at high I is indicative 
of what energy injection might do. The lower i power is sensitive to 
whether nearby clusters are in the map or not, more likely in the 25 
deg^ case than in the 4 deg^ case because of smaller sample variance. 
As well, all large-scale tidal power is included in the peak-patch algo¬ 
rithm, whereas the TreeJ-SPH periodicity results in power truncation 
at the fundamental mode of the box. 
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Figure 4. The catalogue of clusters up to z = 2.5 in a 2° x 2° pyra¬ 
mid pseudo-realization generated from our 200 Mpc ACDM simulation. 
The radial direction is comoving position. In (a), colour and size vi¬ 
sualize cluster/group mass and size found with a halo identification 
algorithm. In (b), colour visualizes the optical detection probabilities 
of the clusters and length visualizes the uncertainty in redshift local¬ 
ization, using parameters appropriate to the 100 deg^ RCS survey of 
Yee and Gladders (2001). Both decrease with mass. For redshifts be¬ 
yond unity the uncertainty in redshift is about the size of the box of 
our hydrodynamical simulations, 200 Mpc, with of course much bet¬ 
ter angular positioning, (c) and (e), the projections of (a) and (b) 
onto 2D maps should be contrasted with the 1.5'-smoothed SZ map 
of the same pseudo-realization: source-blending due to superpositions 
is evident. For clarity, only clusters of M > 10^^'® Mq are shown 
in (a) and (c) and only clusters of non-zero detection probability are 
shown in (b) and (e). In spite of the RCS uncertainties, it is clear that 
such optical surveys can provide powerful templates for source-finding, 
source-deconvolution, and redshift localization. 
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olution hydro simulations with entropy feedback and cooling included to refine 
our SZ predictions. For the adiabatic calculations shown here, we have demon¬ 
strated that, although gas distributed in filaments and membranes in the IGM 
have non-negligible SZ signals, the observable range for planned experiments 
is dominated by the “extended-source” picture. We also addressed the issue 
of source identification in SZ maps of the size, resolution and noise level that 
mimic planned observing campaigns, and explored how optical photometric sur¬ 
veys such as the “Red-Sequence Cluster Survey” can be used to help to find SZ 
sources and localize them in redshift. This could allow us to identify superclus¬ 
ter structures, i.e., to deproject SZ maps. Indeed, given that large area optical 
surveys will be available, an argument can be made for catalogue-targeted rather 
than blank-field SZ observations. We expect both strategies will be used. 

To understand the complex residual maps after unwanted extended-source 
removals, a large number of fast Monte Carlos will be essential, precluding a 
hydro-only approach with current computational power. The extended-source 
approach applying constrained mean pressure-profiles to halos identified in ultra- 
large low-resolution N-body simulations or by the peak-patch technique, cali¬ 
brated with the high resolution hydro results, is a worthwhile path to further 
explore. The snapshot of the work described here mostly shows that there is 
much to develop in SZ numerical simulations and in analysis tools to directly 
confront the exciting new SZ datasets coming down the pipe. 
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